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We show how to efficiently exploit decoherence free subspaces (DFSs), which are immune to 
collective noise, for realizing quantum repeaters with long lived quantum memories. Our setup 
consists of an assembly of simple modules and we show how to implement them in systems of cold, 
neutral atoms in arrays of dipole traps. We develop methods for realizing robust gate operations on 
qubits encoded in a DPS using coUisional interactions between the atoms. We also give a detailed 
analysis of the performance and stability of all required gate operations and emphasize that all 
modules can be realized with current or near future experimental technology. 
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I. INTRODUCTION 

One of the major obstacles in the realization of quan- 
tum information processors (QIP) is decoherence, caused 
by the unwanted interaction of the system with its en- 
vironment. Even the feasibility of simple special pur- 
pose quantum devices, such as quantum repeaters, is sig- 
nificantly affected by the presence of decoherence and 
a special precaution to counteract the resulting loss of 
quantum information has to be taken. A number of con- 
cepts have been developed to fight decoherence and in- 
crease the reliability of QIP including active error correc- 
tion 2ii M ^^'^ fault tolerant quantum computing jy] . 
Aside from this it was shown that, dependent on the type 
of interaction processes between the qubits and their en- 
vironment, subspaces of the system's Hilbert space can 
exist which are immune to decoherence processes induced 
by this interaction 0, IE S 111 ■ These decoherence free 
subspaces are a method of passive error correction or er- 
ror prevention and can significantly increase the lifetime 
of quantum information and reliability of QIP as already 
demonstrated with ion traps , N MR and optical meth- 
ods 9, Tol [H El [H El [H [l3 . However, this robust- 
ness makes states within the decoherence free subspace 
(DFS) difficult to manipulate via controlled operations 
and, in general, complicated interactions with a DFS are 
necessary for realizing universal QIP j^j. 

A particularly promising application of DFSs are quan- 
tum repeaters |l7l llSj which require very long quantum 
information storage times but only a limited set of quan- 
tum manipulations. In this paper we show that the quan- 
tum operations necessary for a quantum repeater can 
be realized on DFS qubits with nearest neighbor hop- 
ping [T9I l2Cl|| and interactions with ancilla qubits pH 
only, thus circumventing the difficulties usually associ- 
ated with the usage of DFS Q . Quantum repeaters are 
used to distribute maximally entangled pairs of qubits 
over long distances which is a necessary prerequisite 
for important applications such as quan tum cryptogra- 
phy |22| and quantum teleportation |23j |. They are thus 
considered one of the most important near future special 
purpose QIP. 

The basic idea of a quantum repeater is to divide the 



transmission line into segments with a length of the or- 
der of the attenuation length of the channel. On each 
segment entangled particle pairs are created and by ap- 
plying entanglement swapping and purification pro- 
tocols 0, I23, HE] entangled pairs of larger distances 
are produced. By repeating these steps a distant pair 
of qubits with high entanglement fidelity is obtained. 
Entanglement purification, distillation, and concentra- 
tion have been realized in a number of seminal experi- 
ments 27, 28, 29, 30] using photonic states for encod- 
ing the qubits. Techniques for converting flying pho- 
tonic qubits into atomic ones via cavity-enhanced inter- 
actions 1^ llll or the usage of atomic ensembles 
have been proposed. Very recently a demonstration of 
long distance entanglement of massive particles has been 
achieved |35l| . 

Here we develop a proposal for a quantum repeater 
with DFS quantum memory using the entanglement pu- 
rification scheme introduced in |l7l l36l | (also known as 
the "Innsbruck protocol" ) . During the purification pro- 
cedure the entanglement of a primary pair of qubits is 
increased by sacrificing auxiliary entangled qubit pairs. 
The entanglement fidelity of the final entangled qubit 
pair does not converge to 1, but to a value which de- 
pends on the fidelity of the initial primary pair and on 
the constant fidelity of the auxiliary pairs. However, the 
scheme 0, Is^l is favorable from a pragmatic point of 
view since it requires fewer physical resources than the 
purification schemes used in [2^ HE 113 ■ On the other 
side, the protocol 0, Is^ takes longer and thus quan- 
tum information has to be stored reliably for the whole 
duration of the protocol. We overcome this apparent 
disadvantage by the usage of DFS qubits of four neutral 
atoms which form the key elements of our DFS quantum 
repeater. 

Our scheme is based on current and near future tech- 
nology and is motivated by the extraordinary progress of 
the field of cold, neutral atoms in dipole traps over the 
past few years both on the theoretical |33, |33, EE 113 ^'^d 
on the experimental |42l l43l | side. Atomic quantum reg- 
isters have been realized in arrays of microlenses ,4^, 
and by exploiting the superfluid to Mott-Insulator transi- 
tion in optical lattices HEEEESj single and many atoms 
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have been manipulated coherently |48l I49| , atoms stored 
in standing waves have been transported over macro- 
scopic distances using conveyor belt techniques |H3, IHH 
Is^. l53l|. and multipartite entangled states of atoms were 



created |49|, |5Jj . Furthermore, extremely versatile tech- 
nologies such as optical tweezers and holographic traps 
that are controlled by liquid crystal spatial light modula- 
tors 55, 56, 57] are currently developed. They allow for 
essentially arbitrary designs of optical potentials. 

We will show that a DFS quantum repeater can be as- 
sembled from modules which only require rotations of the 
logical DFS qubits around the x and z axes by angles tt 
and 7r/2 and a controlled phase (CPHASE) gate where an 
ancilla atom acts as the control qubit. The single qubit 
rotations will be realized by selectively lowering poten- 
tial barriers between atoms and the system will not leave 
the DFS during these gates. The ancilla qubit, which is 
required to interact with a DFS qubit for realizing the 
CPHASE gate, will not be decoherence free. However, 
we will show that its influence on the DFS memory is 
sufficiently small to allow for gate fidelities better than 
98.7%. Furthermore, we will show that our scheme allows 
for deviations of the system parameters of a few percent 
without significantly affecting the gate fidelities. 

This paper is structured as follows: In Sec. ^wc will 
briefly summarize the repeater protocol and show how 
DFSs can be embedded into this scheme. The goal of 
this section is to show what resources are needed for the 
implementation of the quantum repeater. In Sec. Illll we 
propose methods of how the operations of Sec. ^ can be 
physically realized with cold, neutral atoms in arrays of 
dipole traps. The results are summarized in Sec. IIVI 



II. QUANTUM REPEATERS WITH 
DECOHERENCE FREE SUBSPACES 

In the following we will describe the purification and 
repeater protocol of Refs. ^17, ,i36j and extend it by us- 
ing DFS. The basic steps of the purification scheme are 
illustrated in Fig. ^ Only one quantum channel (repre- 
sented by the lines with an arrow) and two DFS qubits 
(represented by the solid circles) are needed. The goal is 
to share a high fidelity Bell state between A and B which 
is done by transmission of "one-half" of a maximally en- 
tangled flying qubit pair (open circles in Fig. ^i)). In 
general, the flying qubit (e.g., a photon) is not directly 
sent from Ato B since the present available transmission 
quantum channels (particularly optical fibers) degrade 
the entanglement of the flying qubit pair exponentially 
with the distance between the nodes. Thus, if the chan- 
nel is too long, the minimum fidelity required for the 
successful application of the purification protocol (i^min) 
is not achieved. To overcome this problem a sufficiently 
high number of intermediate nodes between A and B is 
introduced such that the distance between two nodes is 
of the order of the optical fiber's attenuation length [36l |. 
Entangled flying qubit pairs are prepared at each of the 
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FIG. 1: Illustration of the entanglement purification scheme. 
Pairs of entangled flying qubits (open circles) are created via 
entanglement swapping. These pairs are then used to purify 
an entangled pair of DFS qubits (solid circles; see text). The 
framed steps (iv)-(vi) are repeated until a desired entangle- 
ment fidelity is achieved. 



intermediate nodes and transmitted to the next node. In 
Fig. n this is illustrated for the example of one interme- 
diate node. After this, entanglement swapping at all 
intermediate nodes is performed which yields a partially 
entangled pair of flying qubits between A and B and its 
state is transferred to a DFS qubit pair (step (ii)). Fi- 
nally, as in (i) and (ii), auxiliary entangled qubit pairs 
are generated which are sacrificed in order to increase 
the entanglement of the DFS qubit pair (step (iv)-(vi)). 
Steps (iv)-(vi) are repeated until a desired degree of en- 
tanglement is obtained for the DFS qubit pair. 

However, a straightforward extension of this scheme in 
order to connect arbitrary separated points A and B is 
not possible since the distance between the intermediate 
nodes has to be chosen such that the fidelity of the fi- 
nal pair is not lower than Fmin- This is, in general, not 
possible if the distance between A and B is too large. 
Therefore, a quantum repeater concept by means of a 
nested purification scheme was proposed |l7l l36j| which 
uses the above purification protocol to create a number 
of entangled DFS qubit pairs which are then used to pu- 
rify a lower number of additional DFS pairs with greater 
distance. These pairs are then again used to purify an 
even lower number of pairs, etc. until merely one entan- 
gled pair is left. This procedure is illustrated in Fig. El 
The boxes stretching over lines (i) and (ii) represent the 
purification protocol of Fig. Q] The entangled pairs of 
line (ii) are then used to purify the pairs in line (iii) in 
the same way as in Fig. ^ except that the flying qubits 
have to be replaced by DFS qubits, i.e., entanglement 
swapping between DFS qubits, a state transfer between 
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FIG. 2: The quantum repeater scheme. The dashed boxes in- 
dicate the subroutine shown in Fig. Q The created entangled 
DFS qubit pairs are used to purify further DFS qubit pairs 
(in line (iii)). The resulting pairs are then used to purify the 
pairs of line (iv) and so on. 



DFS qubits (to initialize the pairs in line (iii) and (iv)) 
and purification by using entangled DFS pairs have to 
be performed. The diagram shown in Fig. El can be ex- 
tended in horizontal and vertical directions such that an 
arbitrary distance between the final pair can be achieved. 
It was shown that the number of (DFS) qubits necessary 
for the scheme increases merely logarithmically with the 
distance of the final entangled pair [itI Is^ . 

In summary, the purification procedure shown in Fig.^ 
can be realized by assembling the following modules: (0) 
The creation of pairs of fiying qubits and performing en- 
tanglement swapping between them, (1) the realization of 
two DFS qubits, (2) the transfer of the state of the flying 
qubits to the DFS qubits, and (3) the actual purification 
of the entangled DFS qubit pair. Additionally, for the 
nested purification protocol, i.e., for the full implemen- 
tation of the quantum repeater, we have to (4) perform 
entanglement swapping between neighboring DFS qubits, 
(5) we need a state transfer between DFS qubits (analo- 
gously to module (2)), and (6) we must perform entan- 
glement purification using two DFS qubit pairs. Module 
(0) can be implemented using standard quantum optical 
methods as demonstrated in 58, 59, 60, 61, 62J. It docs 
not involve DFS qubits and will thus not be discussed 
further in this paper. 

In the following we will show how each of the remaining 
modules can be realized. For the operations using ancilla 
qubits and thus leaving the DFS we then calculate the 
achievable fidelities. 



A. Realization of the repeater modules 

1. Module 1: Decoherence free memory and single qubit 
rotations 

The main sources of decoherence in systems of cold 
atoms are fluctuations of the opti cal p otential and exter- 
nal electromagnetic (stray-)fields [63. These sources are 
homogeneous across the region of the physical qubits and 
the interaction with them thus only depends on their in- 
ternal states. The Hamiltonian describing the situation 



Here, i labels the physical qubits, B = {B^, By, B^) de- 
scribes the fields and fiuctuations, and the cr^ are Pauli 
matrices acting on the ith physical qubit. 

The corresponding decoherence free subspace is 
spanned by the states with total zero angular momen- 
tum, and the lowest number of physical qubits giving a 
twofold degenerate DFS is 4. The DFS qubit is then 
represented by 8, 64] 

|0)i = J(|01)-|10))®(|01)-|10)), 

(2) 

ii)^ = ^(2 11100) -t- 2 10011) - (101) + mr') , 

where \ijkl) = K)i |j)2 1^)3 IO4 with i,j,k,l = 0,1 de- 
scribes the computational basis states of the four physical 
qubits. Each of the four physical qubits is represented by 
two internal states of an atom which are trapped in an 
array of dipole traps, for example, an optical lattice. By 
the method described in Sec. lIII j^ the system is prepared 
in the logical state |0)^ and in Sec. lIIIBl it is shown that 
by tuning the system parameters appropriately we can 
implement the Hamiltonians Hx and Hz with 

Hx = Xl = ■^iVi2 + 2V23), (3) 

Hz = Zl = -Vi2, (4) 

where Xl {Zl) is the X-gate (Z-gate) for the logical 
qubit and Vij are permutations of the physical qubits i 
and j. Note that the logical Z-gate can also be repre- 
sented as Zl = —V^i. Since the DFS defined by Eq. 
is invariant under permutations of the physical qubits j3| , 
the time evolution of the system generates rotations 

RliO) - e-i^-, (5) 
Rl{e) = e-'i^^ (6) 

without ever leaving the DFS, where 9 = 2t and t is the 
time necessary for performing the rotations. The index L 
indicates that these are operations on the logical qubit (in 
the following we will also use this notation for operations 
on physical qubits, where we omit the index L). 

Although it is not required in the repeater protocol, we 
can thus perform arbitrary rotations of the logical qubit 
(and hence we can prepare arbitrary states in the DFS) 
since every two level unitary operator can be decomposed 
in Rl and Rl 

2. Module 2: State transfer between flymg and DFS qubit 

The quantum circuit for the state transfer of the fiying 
qubit to the logical qubit is shown in Fig. O The state 
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|0)^ - Rl{^) - -Z^ - Rl{-%) 



FIG. 3; Quantum circuit for the state transfer of the auxihary 
qubit to the DFS qubit. 



of the flying qubit is denoted as \ip)^ — a |0)^ + (3 
As described in Sec. IIIII the state of the flying qubit 
(a photon) is first transferred to an appropriate inter- 
nal state of an auxihary atom. The state of this atom 
then controls the interaction between physical qubits 3 
and 4 leading to a controlled SWAP operation between 
these atoms. Since = —V34, a controlled SWAP oper- 



ation of the physical qubits 3 and 4 leads to a controlled 
{—Zl) (CPHASE) gate on the logical qubit. As can be 
seen from Fig. |31 apart from rotations given in Eqs. © 
and we need to perform a rotation R^{—n/2) and 
a measurement of the control atom which can be done 
by standard techniques such as laser induced Rabi oscil- 
lations and state selective laser excitations. Eventually, 
the state of the atom is transferred to the logical qubit, 
i.e., |V)i=a|0)^+/3|l)^. 



3. Module 3: Entanglement purification 

The local operations required for entanglement purifi- 
cation are shown in Fig. 0] 



auxihary qubit — U^'^ — R^i 







Rl{-^) 





FIG. 4: Quantum circuit for entanglement purification. The procedure is performed on both sides, A and B. 



This circuit corresponds to the scheme proposed in |2fi| . 
but here it is expressed in terms of operations that we can 
implement in our setup. The circuit has to be applied on 
both ends of the transmission line A and B, where 



(7) 



As in subsection III A 21 the auxiliary qubit is again an 
atom that controls a phase operation on the logical qubit. 
If the outcome of the measurements at A and B coincide 
we end up with an entangled DFS qubit pair of higher 
fidelity. If they do not coincide the purification procedure 
has to be repeated. 



4- Module 4: Entanglement swapping between DFS qubits 



Entanglement swapping between neighboring DFS 
qubit pairs can be done by a Bell measurement and sub- 
sequent local operations i?|(7r/2) and i?|^(7r/2) that are 
controlled by the outcome of the measurement l3^ . The 
Bell measurement can be decomposed in a controlled not 
(CNOT) gate, a Hadamard gate, and a measurement in 
the computational basis (the Hadamard gate can be re- 
placed by, e.g., the sequence i?£(7r/2)i?£(7r/2)i?£(7r/2)). 
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FIG. 5: Quantum circuit for a CNOT operation between two logical qubits. 



The implementation of a decoherence free CNOT opera- of two-particle interactions with pairwise individual tun- 
tion is relatively complicated including lengthy sequences 
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Therefore we use an ancilla atom (as in Sec. lII A'2l 
and 111 A'S)| to mediate the interaction between different 
DFS qubits. A possible quantum circuit for a CNOT gate 
between DFS qubits is shown in Fig. [S] The centerline in 
the circuit represents the ancilla qubit. The shown cir- 
cuit uses only gates that we can implement in our scheme 
and leads to a CNOT gate with the upper DFS qubit as 
the control qubit. We will discuss in Sec. Ill Bl how the 
decoherence of the auxiliary atom influences the fidelity 
of this operation. 

For the measurement of the DFS qubit we use the same 
ancilla as for the CNOT gate. The corresponding circuit 
is shown in Fig. |51 The final measurement of the ancilla 



11), ^i^-(f) 



-Zl 



FIG. 6: Measurement of the logical qubit. 
is equivalent to a measurement of the DFS qubit. 



5. Modules 5 and 6: State transfer between DFS qubits and 
entanglement purification using pairs of DFS qubits 

The state transfer of one DFS qubit to another can be 
done analogously to Fig.O where has to be replaced 
by . The necessary CPHASE gate between DFS 

qubits can be constructed similarly as the CNOT gate: 
We merely have to omit the single qubit rotations in the 
last line of Fig. Also the entanglement purification 
scheme of Fig. 0] can be adapted in this way. Again, 
the auxiliary qubit is replaced by a DFS qubit and the 
CPHASE (or CNOT) gate is realized as described above. 



B. Operation time and decoherence 

Some of the above operations require the use of a single 
auxiliary qubit which is not decoherence free. Entangling 
this qubit with the DFS qubits will therefore lower the 
overall fidelity of the corresponding gate operations on 
the DFS qubits. The auxiliary qubit is represented in 
our scheme by two long-lived internal states of an atom 
and thus the major experimentally expected decoherence 
mechanism is dephasing |66j |. This process is described 
by the master equation 



dtQaux — 



^ < 

2 



c) , (8) 



where ^?aux is the density matrix of the auxiliary atom, H 
is the Hamiltonian of the system, I/7 characterizes the 
decoherence time, and is a Pauli matrix. Decoherence 
of the DFS qubits is not taken into account here. 



A quantitative measure of the operation fidelity is de- 
rived from Ref. A distance measure of two density 
matrices a is given by 



!{b, ct) = TrA/ 



(9) 



where, in contrast to |65|, we use the squared trace since 
for pure states this corresponds to a probability, i.e. the 
absolute square of the overlap of the two states. The 
operation- (or gate-) fidelity can be derived from this 
quantity: Let f be a quantum operation describing a 
quantum circuit where decoherence is taken into account, 
and £q the same operation without decoherence. We de- 
fine the operation fidelity by 



F(£o,f) =min/(£o(g),f(f?)) 

g 



(10) 



where it is sufficient to minimize over pure input 
states Q 

Let us first consider the state transfer from the atom 
to the DFS qubit (Fig.|3Jl. To minimize decoherence ef- 
fects, the i?£(7r/2) and Rf^{—n/2) gates are appfied be- 
fore the auxiliary atom interacts with the DFS qubit and 
after the measurement, respectively. In addition, we ne- 
glect decoherence during the i?^(— 7r/2) rotation since the 
time to perform such a gate on one atom is, in general, 
much smaller than the expected decoherence time I/7 
(see Sec, nil D|) . Thus, decoherence occurs mainly during 
the CPHASE gate. We examined the dynamics of the 
system during the application of the CPHASE gate by 
solving Eq. © analytically. The resulting gate fidelity 
for the CPHASE operation is given by 



-Fcphase(^-o,'^^) 



1 



-7t 



> 1 



2 



(11) 



As will be discussed in Sec. lIII Dl this behavior is sufficient 
to ensure a working repeater scheme. 

Note that, since unitary operations do not change the 
fidelity of the circuits [s^, all operations which only in- 
volve one CPHASE gate have the same fidelity given in 
Eq. (|ll|l . This includes particularly the operations shown 
in Figs. Hand El 

In a similar way we calculated the fidelity of the CNOT 
gate between the two DFS qubits (Fig. |SJl. Again, all 

i?^^'^'' gates are applied before the auxiliary atom inter- 
acts with the DFS qubit and after the measurement, re- 
spectively. During these periods the auxiliary atom is not 
entangled with either DFS qubit and thus its decoherence 
has no infiuence on them. Furthermore, we assume again 
that decoherence during the rotations can be ne- 

glected. The gate fidelity of the CNOT operation is then 
given by 



^CNOt(^-0!'^^) 



1 + c 



-'ft 



(12) 



where t is the time which is needed to perform one 
CPHASE gate. Since the two DFS qubits are not en- 
tangled before the CNOT gate is applied it is sufficient 



to use only pure, unentangled states for the minimiza- 
tion [ill. The fidelity for a CPHASE gate between two 
DFS qubits is the same as for the CNOT gate because 
they differ only by unitary operations. As in the case 
of the CPHASE gate, the decay of the CNOT fidelity is 
slow enough to ensure a working quantum repeater (see 
Sec lmPll . 

The fidelity of the state transfer between two DFS 
qubits described in Sec. Ill A 5l can also be calculated by 
explicitly solving the time evolution. The state transfer 
fidelity Fst is given by 

Fst{£o,£) - J (1 + c-'^*) (1 + c-"*) ^ 1 - • (13) 

The fidelity of the entanglement purification using 
pairs of DFS qubits (see Sec. Ill A 5() is calculated using 
the chaining property for fidelity measures |65l | . Assum- 
ing that the two DFS qubits on one side A oi B are not 
entangled with each other we find as a lower bound 

FEp(fo,f) > 1- Q + V2^7i- (14) 

In both formulas, Eqs. H13() and H14() . the time t is again 
the time needed for performing one CPHASE gate. As 
we will show in Sec 1111 D1 all of these infidelities are small 
for typical experimental decoherence times. 

In summary, we have shown in this section that the 
DFS quantum repeater can be divided into simple mod- 
ules which require (apart from single qubit rotations and 
measurements of an ancilla atom which can be realized 
by standard methods) merely the operations i?2(±^), 
i?£(±7r), i?|,(±f), RUtt), and a CPHASE gate on a 
DFS qubit with the ancilla atom as control qubit. In the 
following section we proceed by describing how these ba- 
sic operations can be implemented in systems of trapped 
neutral atoms. 



III. PHYSICAL IMPLEMENTATION IN 
ARRAYS OF DIPOLE TRAPS 

The physical system we consider for the implementa- 
tion of a DFS qubit is a one-dimensional chain of four 
atoms in an array of optical dipole traps, e.g., an optical 
lattice with one atom per site. Each of the atoms 
has two long-lived internal states a and b which encode a 
physical qubit |vac) = |0), w |vac) = |1). In Fig.[7|two 
DFS qubits are schematically shown. The thick, diago- 
nal line represents an optical conveyor belt 53] storing 
the auxiliary atom. This atom can be transported to 
the DFS strings where it occupies an empty site between 
the third and the fourth atom, which is necessary for 
the CPHASE operations described in Sec.|H] Aside from 
this, the auxiliary atom is needed for the state transfer 
of the flying qubit to the DFS qubit: The flying qubits in 
our repeater scheme are represented by entangled photon 
pairs. After a photon arrives at a repeater node its state 
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/ 

FIG. 7: Schematic picture of a repeater node. The four 
register atoms per qubit (black circles) are stored in a one- 
dimensional array of dipole traps, e.g., an optical lattice. Be- 
tween the third and fourth atom there is a free lattice site in 
which the auxiliary atom (gray) can be moved, e.g., via an 
optical conveyor belt. 

is transferred to the auxiliary atom, e.g., by a cavity- 
enhanced interaction |^ [s^ • The atom is then moved 
to the DFS qubit and the operation shown in Fig. [3| is 
performed. In order to implement the rotations on the 
DFS qubits we also have to be able to selectively lower 
lattice barriers between certain lattice sites. This can be 
done by using suitable superlattices or by applying sin- 
gle laser pulses as described in The addressability 
of single atoms is only required for the initialization of 
the register, cf. Sec. IllI Al 

The general Hamiltonian for a one-dimensional atom 
chain for M lattice sites (without auxiliary atom) is given 

by m 

(15) 

where jj"-* are state dependent hopping matrix elements 
from site j to site j + \ and f/a,b, Uab are state dependent 
on-site interaction strengths. We should mention that for 
the rotations i?£ and and for the CPHASE gate the 
hopping constants for both kinds of atoms are required 
to be the same, = Jj''\ as well as Ua = Ub = Uab, 
cf. also Appendix|Al Apart from the initialization of the 
register we need M — 5 lattice sites. 

In the following subsections we will show how the reg- 
ister can be prepared in state |0)^ and how the rotations 

r'i'''^ and the CPHASE gate can be implemented. These 
three ingredients are sufficient to realize the repeater pro- 
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tocol. 



1. Analytical treatment 



A. Initialization of state |0)^ 

The preparation of state |0)^ requires the creation of 
a pair of sing let states (|01) - |10))/V2. These states are 
decoupled from global operations but can be realized by 
local operations and by selectively lowering the potential 
barriers between neighboring atoms. Starting out from 
the product state |0110) the potential barrier between the 
first and second atom and between the third and fourth 
atom are lowered (e.g. via a superlattice potential). This 
generates two double well type potentials. The parame- 
ters of the Hamiltonian Eq. (|15|l are then given by 



^2=0. 



"^1,3 



T(f') 

"^1,3 



J, Ua^Ub^ Uab 



U. 



(16) 



Assuming that J/U is sufhciently small (see also 
states with two atoms in one site can be adi- 



Sec. lmBt 

abatically eliminated leading to an effective Hamiltonian 
El of the form 



rV^y 



alal + alal). 



(17) 



Turning this Hamiltonian on for a time t = irhU/SJ'^ we 
obtain the state (|01) + i |10)) (|10) + i |01))/2. Finally, 
a local operation on the first and fourth physical qubit of 
the form exp(— iTrtrf /4) and exp(— i7r(7|/4), respectively, 
yields the desired product of singlet states. Alternatively, 
methods proposed in can be used to realize the state 
(|01) + |10))(g) (|10) + 101))/2. Local single qubit rotations 
of the first and fourth qubit can be used to bring this 
state into the desired from. If laser addressing of single 
physical qubits in adjacent lattice sites is not possible 
we can move the first and the fourth atom into their 
neighboring, free sites (and back again) by lowering the 
corresponding potential barriers. 

In order to write information into the register we re- 
quire a free lattice site between the third and fourth 
atom. Hence, in a last step the fourth atom is moved 
into the fifth lattice site by lowering the potential barrier 
between the fourth and fifth site. 



B. Rotations of the logical qubit 

In this section we analytically investigate how to imple- 
ment the rotations i?£ and Rf^ in the two limiting cases 
of a very high interaction between the register atoms and 
a vanishing interaction. The system is simulated numer- 
ically to check our results. 



The Hamiltonian Hz is, apart from a sign, equivalent 
to a swap of the first two atoms V12 [see Eq. Q]. If we 
assume high potential barriers between the second, the 



third, and the fourth atom, i.e., J2 



(a,6) 



J. 



(a,6) 



0, such 



that these atoms form a separate, decoupled system, we 
merely have to describe the dynamics of the first and 
the second atom. If we further assume that Ua,b, Uab are 

much larger than the hopping terms jj"'''^ we can adia- 
batically eliminate all states with more than one atom in 
a lattice site. By choosing the parameters according to 



4''^ = jf ' = J, 

Ua = Ub = Uab = U, 

we arrive at a the effective Hamiltonian 

Hz = -— {(Ji(J2 + + ^1^2 + I12) 



2J2 
U 



2J2 

17" 



(18) 
(19) 



(20) 



where we omitted a constant term. The time evolution 
of the system then generates the rotation 



Rl{Q) = exp 



(21) 



where d — A,Pt/Uh. Hence, the time for a i?£(7r)- 
rotation is given hy t = TrUh/4J^. As an example we 
consider Na atoms in an optical lattice with a wave- 
length of A = 514nm. The height of the transverse po- 



tential barriers are assumed to be V,, 



which is in the range of today's experiments [43|. Here, 
= h?{2Tr)^ /2mX^ is the recoil energy and m is the 
mass of the atoms. If we increase the s-wave scatter- 
ing length by a factor of 5 using a Feshbach resonance 
and choose Vx = 7.7E^ between the first and second lat- 
tice site we get a ratio U/J — 75 which is well in the 
regime where the adiabatic elimination is valid, as will 
be discussed in more detail below. Note that for these 
parameters the assumptions for the validity of the Bose- 
Hubbard Hamiltonian Eq. (|15|l are fulfilled. The hopping 
term is then calculated to be J = 0.033i?R which results 
in an operation time oit = 8.7ms for the i?|^(7r)-rotation. 
We emphasize that we never leave the DFS during the 
rotation and it is expected that this time is much shorter 
than the decoherence time of the DFS qubit. 

Nevertheless, the time for a i?£(7r)-rotation can be 
made considerably shorter by tuning the interaction be- 
tween the register atoms to zero (e.g., via a Feshbach 
resonance). As we show in Appendix IbI the qubits still 
remain in the DFS for all times. In this interaction- 
free case, the exact dynamics of the two atom system 
(i.e., without adiabatic elimination) results after a time 
t = TTh/2J in the mapping 



|a>/3)i2 



1/3, a) 12 , a,/3 = 0, 1 



(22) 
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FIG. 8; Numerical studies of the phase difference (j> between 
states |0)^ and after applying the rotation Rtiji) versus 
the ratio U/J calculated with the full Hamiltonian Eq. 
and the parameters Eqs. 1181 and 1191 given by adiabatic 
elimination. 



Up to a global phase this is equivalent to a i?|^(7r)- 
rotation on the logical qubits. The time needed for a 
lattice as above \s t = 0.23ms. The disadvantage of this 
scheme is that only rotations with a phase of tt can be 
performed. 

According to Eq. (O the implementation of Hx in- 
volves three atoms. In the same spirit as for the Z^-gaXe 
we require that the potential barrier between the third 



and fourth lattice site is high, i.e. 
ing the remaining parameters as 



J. 



(a,b) 



0. By choos- 



J^"' = \/2 j}'^^ = V2 ) = = V2 J , (23) 

Ua = Ub = Uab = U , (24) 

and requiring that J /U 1, we again adiabatically elim- 
inate all states with more than one atom in one well. A 
detailed calculation is given in Appendix^ The result- 
ing Hamiltonian on the logical qubits is then given by 



2V3J^ 
U 



Xr 



(25) 



where we omitted a constant term. This Hamiltonian 
generates a rotation 



i?|(e) = exp ( -i-Xi 



(26) 



where 9 = —i^/SJ^t/UH. Taking the numbers of the 
above example, the time required for a i?£(— tt) rotation 
is i = 5.0ms. 



2. Numerical simulation 



To test the validity of the adiabatic elimination we sim- 
ulated the system numerically. In case of _R|^-rotations 
the parameters of the exact Hamiltonian Eq. (|15|l are 




FIG. 9: Fidelity Eq. ^ for the rotation RUtt). In (a) the 
hopping term between the second and third atom is detuned, 
in (b) the interaction between register atoms of kind a and 
b. The other parameters are (a) Ua,b = Uab = 100 Ji, (b) 

Ua,b = U ^ WOJ. 
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FIG. 10: Numerical studies for a fast rotation with no or only 
weak interaction between the atoms. The fidelity Eq. H27|l 
of the rotation -R|,(7r) and the phase difference between the 
two logical qubits after this rotation versus the ratio [// J is 
calculated with the full Hamiltonian Eq. 1151 . The accuracy 
of the gate gets better with a smaller interaction U between 
the atoms. Note that here no adiabatic elimination has been 
used. 



chosen according to Eqs. H18() - (|19() .We numerically cal- 
culated the fidelity of the DFS qubit remaining in its 
state during an _R£-rotation, 



/z^ |(Z|^exp(-iiJztA)|/)j2 



(27) 



with / = 0, 1. For U/J — 75 it turns out that the in- 
fidelity 1 — is smaller than 3x10^'^ for / = 1 and 
smaller than 10^^ for I = during the whole duration of 
a i?£(27r)-rotation. 

Since this fidelity does not take the phase difference 
between states |0)^ and |1)^ into account, was calcu- 
lated as well. In Fig. |S1 the phase difference is shown for 
the rotation i?|^(7r). As expected, for increasing ratios 
U/J the phase difference gets closer to tt. For U / J ~ 75 
the inaccuracy of the achieved phase is already smaller 
than 10~^, which shows that the adiabatic elimination is 
well justified. 

From an experimental point of view it is also impor- 
tant to know what happens when some of the parameters 
are detuned. Some results are shown in Fig. |51 The error 
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FIG. 11: Numerical studies of the fidelity Eq. ^ of the 
rotation RKn) versus the ratio U/J calculated with the full 
Hamiltonian Eq. and the parameters Eqs. <\2'A^ and H24|l 
given by adiabatic elimination. 




FIG. 12: Fidelity Eq. (EHl for the rotation Rl{n). In (a) the 
hopping term between the second and third atom is detuned, 
in (b) the interaction between register atoms of kind a and 
b. The other parameters are (a) Ua,b = Uat = 100 Ji, (b) 

Ua,b = U = WOJ. 



due to the nonvanishing hopping between the second and 
third lattice site J2 is negUgible as long as J2 < 10^^ Ji. 
Because of this we have completely neglected the influ- 
ence of the fourth register atom which is not involved in 
these rotations. In the case of a detuning of the inter- 
action coefficient Uab the fidelity changes quite quickly if 
Uab/U differs from 1, but is smaller than 10"'^ as long as 
the detuning is less than 2%. 

The fidelity Eq. I|27|) and the phase difference ip be- 
tween the logical qubit for the faster scheme of a i?£(7r)- 
rotation for small, residual interactions U /J (in the ideal 
case this interaction is vanishing) are shown in Fig. 1101 
For this case the inaccuracy of the gate decreases with 
decreasing U/J and is smaller than 2 x 10"'^ for U/J w 
10~^. The error calculated according to Eq. H27|) is orders 
of magnitudes smaller than the phase inaccuracy. 

The simulations of the i?£-rotations lead to similar 
numbers. The fidelity of this process, determined by 

/x^|(0|^exp(-ii/xtA)|l)iP, (28) 

is shown in Fig. [TTl depending on the ratio U/J. As can 
be seen from this figure, the infidelity 1 ~ fx decreases 
with increasing ratio oi U/J and for a [// J w 75 the 
infidelity is smaller than 10"'^. Any infidelities because 
of a change of the relative phases between the two logical 
states are orders of magnitudes smaller and play no role. 

Figure IT^ shows the results of numerical simulations if 
various system parameters are detuned. For the case of 
a nonideal hopping the deviation should be smaller than 
1.5% to yield infidelities smaller than 3 x 10"'^. A slight 
detuning in the interaction strength is less critical, but 
should nevertheless be smaller than 4% to achieve small 
errors. 



C. Controlled phase gate 

The CPHASE gate can be realized by a controlled swap 
operation between the third and the fourth atom since. 



as mentioned in Sec. Ill A II = — V34. We assume 
that the potential barriers between the first, second, and 
third atom of the DFS qubit are high enough to ignore 
the dynamics of the first two atoms. The auxiliary atom 
is assumed to be distinguishable from the DFS atoms 
and is confined in a sufficiently deep lattice such that it 
always remains at its site. If the auxiliary atom is moved 
to the free site of the DFS string (cf. Fig. [TJ and if a 
register atom tunnels into the site of the auxiliary they 
interact with each other giving the Hamiltonian 



^CPHASE = H+ Ulala^qUa + UfblhqUa , (29) 

cr=0,l 



where H is the Hamiltonian given by Eq. (|15(l . is 
the state-dependent annihilation operator of the auxil- 
iary atom with q its internal state, and U^^ the interac- 
tion strength which is independent of the internal state 
of the DFS atoms. If the interaction Uq for state |0)^ of 
the auxiliary atom is very large compared to the hopping 
strength, the third or fourth register atoms cannot en- 
ter the lattice site of the auxiliary atom and hence they 
will remain in their sites. Numerical tests show that dur- 
ing the gate operation this still holds for an interaction 
strength of C/q/J ~ 100, where J is the hopping term of 
the register atoms. Under these circumstances the prob- 
ability of finding the register in its original state is larger 
than 1 — 10"'^. If the auxiliary atom is in the state |1)^ 
and Uf — the existence of the auxiliary atom can be 
ignored and we only need to solve the dynamics of two 
register atoms in three lattice sites. A similar scheme has 
been discussed in [6^ [t^ , where an electromagnetically 
induced transparency-like configuration was used to con- 
trol the interaction strength. We will investigate in the 
following the two limiting cases of very large and vanish- 
ing interaction strength between the register atoms. 
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FIG. 13: (a) Numerical results for the fidelities Eqs. (131)11 
and Ij31|l . where the auxiliary atom has a nonvanishing rest 
interaction with the register atoms, (b) Numerical results for 
the fidelity Eq. I|35|l due to the phase difference, where the 
auxiliary atom has a nonvanishing rest interaction with the 
register atoms. In both cases we took for the register atoms 
U/J = 100. 

1. Large interaction between the register atoms 

For the case of large U/J we calculate the fidelities of 
the swap operation, 

foi = \{0<jl\U{t)\laO)\\ (30) 
foo^ fii^\{OaO\U{t)\OaO)\^ (31) 

numerically, where t is the numerically obtained time for 
the CPHASE gate and U{t) = exp(-iiJcPHASEV^) is 
the time evolution operator corresponding to Hamilto- 
nian Eq. 129|) . The states \aaP) , a,/3 = 0,1 denote 
the states of the third and the fourth DPS atom, \a) 
and I/?), and \a = 1) is the internal state of the auxiliary 
atom, located in between the two DFS atoms. We choose 
J(°'^) = J and Ua,t = Uat = U. 

In Fig.lT^ the fidelities Eqs. (|^ and are shown 
for a small but nonvanishing interaction Ul between the 
auxiliary atom and the register atoms for a ratio U/J = 
100. As can be seen the infidelities are smaller than 10"'^ 
as long as the interaction fulfills C/'/J < 0.05. A further 
numerical observation is that the value of /oo (/oi) for 
very small U^/J gets smaller with increasing U/J. 

The fidelities Eqs. (|30|l and (|31|l do not take into ac- 
count the occurrence of phases. In general, the dynamics 
of the system will yield a phase, 

|acr/3) ^ e'^°'' |/3cra) (32) 

with a, /? = 0, 1. Since the interactions t/J do not depend 
on the state of the DFS atoms we have ipoi = ipio and 
foo = fii and a nonvanishing U^ leads to the behavior 

ZY(t)|0)^=e'^-|0)^ (33) 

Z^(t)|l)^=e'^-^[2(|1100) + 10011)) 

-e'^(|01) + |10))®2], (34) 

of the four DFS atoms, i.e., to the presence of a nonzero 
phase if = ipoQ — ifiQ. Thus the state U{t) |1)^ is not 



FIG. 14: (a) Fidelity after evolving the state for a time t = 
TV h/V2J for the case of noninteracting register atoms, (b) 
Fidelity Eq. (1351 of state |1)^ after the time evolution, where 
a perfect fidelity of the atom swap operation has been assumed 
and the error is only due to the phase difference, see text. 

entirely contained in the DFS anymore. To study this 
error we define the fidelity 

u.-\{i\Lm\i),\'-'-^^^, (35) 

which is shown in Fig. 113b depending on Uf / J . The infi- 
delity 1 — /ph is typically one order of magnitude smaller 
than the infidelities shown in Fig. 113b . 

The time which is needed to perform the swap gate 
has been calculated numerically. For a ratio oi U/J = 75 
and the same lattice parameters as for the i?£-rotations 
we get t = 11.4ms. 

Numerical tests show that for the same ratio U/J — 75 
the probability of finding the qubit outside the DFS is 
smaller than 5 x 10~^ during the whole gate operation 
time. 



2. No interaction between DFS atoms 

A substantially shorter gate time can be achieved for 
the case of a vanishing interaction between the DFS 
atoms. In Appendix ^ it is shown that the qubits still 
remain in the DFS during the gate operation. By ana- 
lytically solving the time evolution of the Hamiltonian it 
can be shown (see also Appendix that after a time 
t = nh/ y/2J the dynamics yields the mapping 

\aa(3) ^\f3aa) , a, 0,1, a =1 (36) 

which is a swap operation between the third and the 
fourth DFS atom and thus a CPHASE gate. For the 
same lattice potential as earlier with Vx — 7.7£'r the 
gate time is given by t = 0.33ms. 

The results of numerical tests for a small rest inter- 
action between the register atoms are shown in Fig. 1141 
Figure [Tlk shows the infidelities derived form Eqs. (|30() 
and (|31|) and Fig. 114b shows the infidelity derived from 
Eq. H35|l all depending on U/J. The fidelities are better 
than 1 — 10^"^ for interactions U/J < 0.05. The results 
of numerical simulations for a nonvanishing rest inter- 
action between the DFS atoms and a state |1)^ of the 
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auxiliary atom yield an infidelity smaller that 3 x 10 
for Ul/J < 0.05. 

In conclusion, a reliable CPHASE gate can be imple- 
mented in both cases. If the auxiliary atom is in the 
state |1)^ and thus a swap of register atoms 3 and 4 
should occur, the interaction [/' between the auxiliary 
atom and the register atoms is required to be smaller 
than 0.05 J. If the auxiliary atom is in the state |0)^ no 
swap should occur. This holds during the gate interac- 
tion time for an interaction Uq larger than approximately 
lOOJ. The difference of more than three orders of mag- 
nitude in the interaction can be achieved either by using 
a Feshbach resonance or by using an electromagnetically 
induced transparency scheme as in [6^. 



D. Gate times and decoherence 

Recent experiments have demonstrated that the deco- 
herence time of quantum information stored in DFSs can 
be increased by several orders of magnitudes compared 
to unprotected qubits. In the case of ion traps decoher- 
ence times of more than 7s have been reported 0| and 
although the DPS scheme in this experiment was differ- 
ent from ours we expect a similar improvement in the 
reliability of our quantum memory. The required storage 
time is found by adding the single qubit and CPHASE 
gate times calculated in previous sections. According to 
these the state transfer from the auxiliary atom to the 
DPS takes 11ms, neglecting the time which is needed to 
implement the rotation and measurement of the auxiliary 
atom. The entanglement purification module described 
in Fig.0]takes 26ms, the CNOT operation of Fig. [Stakes 
28ms. The necessary readout process of the DPS qubits 
takes 11ms. In comparison to these times the additional 
operations in steps (i), (ii), and (iii) of Pig. ^ are ex- 
pected to contribute a i iegligib le amount of extra time to 
the overall scheme IS El El ES El i3 ■ AU these 
times are very short compared to the expected decoher- 
ence time of the DPS qubit. Since the system remains 
in the DPS during the entire single qubit operations we 
conclude that decoherence of the DPS qubits plays only 
a minor role. 

However, in the modules which involve a CPHASE 
gate between the auxiliary qubit and a DPS qubit the 
decoherence of the auxiliary atom has to be taken into 
account since, in general, it gets entangled with the DPS 
qubit thus lowering the overall fidelity. As discussed in 
Sec. IiTbI the fidelity for, e.g., the CNOT gate is better 
than 1 — "ft, where 7 is the inverse decoherence time and 
t is the time to perform a CPHASE gate. In our fastest 
example (without interaction between register atoms) a 
CPHASE gate takes approximately 0.33ms. In Ref. |^ 
decoherence times of 146ms for ~ 50 Cs atoms stored 
in a red-detuned standing wave have been reported. In 
the experiment presented in Ref. E3 the influence of a 
moving conveyor belt is investigated. It is found that the 
transport procedure decreases the coherence time by ap- 



proximately a factor of 2. This would still be long enough 
for our purposes. Por the entanglement purification mod- 
ule, which has the worst fidelity of all our models, the 
influence of the auxiliary atom leads to a fidelity of bet- 
ter than 98.7%. Por the other modules the fidelity due 
to decoherence of the auxiliary atom is better than 99%. 
The usage of a blue-detuned trap could give further con- 
siderable increase of the decoherence time. Aside from 
this, it should be emphasized that all these times have 
been measured for many atoms in a trap. In the case of 
single atoms longer decoherence times can be expected. 
Thus, the influence of the auxiliary atom on the modules 
described in Sec. Ill Al is small and only slightly decreases 
their fidelity. 



IV. CONCLUSION 

In this paper we developed a scheme for the robust 
implementation of a quantum repeater. The qubits at 
the repeater nodes are encoded in a DPS which increases 
the decoherence time of the quantum memory consider- 
ably and thus improves the reliability of the scheme. We 
showed that the quantum repeater can be divided into 
simple modules whose implementation only requires a 
small number of different quantum gate operations. The 
single qubit gates can be realized by selectively lowering 
potential barriers between the atoms without ever leaving 
the DPS. Controlled operations are mediated by a single 
auxiliary qubit which is not decoherence free. However, 
the influence of this qubit on the gate fidelities was thor- 
oughly investigated and shown to be sufRciently small 
as not to significantly affect the quality of the repeater 
scheme. 

Furthermore, we presented proposals for implementing 
the DPS qubits and all necessary operations in systems 
of neutral atoms trapped in arrays of optical dipole traps. 
The times for implementing these gates were calculated 
and shown to be much smaller than the expected deco- 
herence time of the DPS qubits. Numerical examinations 
have shown that the achieved gate fidelities can be bet- 
ter than 98.7% and that the gate operations are stable 
against possible, small deviations of the system parame- 
ters. 

Our scheme is extendable to general purpose quan- 
tum computing with DPS qubits, since the elementary 
one-qubit rotations i?£ and a two-qubit gate, which 
together form a universal set of gates, can be realized. 
However, during the CPHASE gate the DPS is left. This 
effect is very small for our quantum repeater scheme and 
the influence of leaving the DPS in other applications 
as well as the scalability of the scheme is the subject of 
further investigations. 

Recent experiments have demonstrated that it is pos- 
sible to implement and control quantum registers with 
a few atoms in optical lattices |7lj or dipole trap arrays 
|45| and that the decoherence times of the auxiliary atom 
are long enough to ensure a successful operation of our 
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scheme [5j, |63| . The experimental combination of atom 
registers and optical conveyor belts is also planned [7l| . 
Hence we are confident that the presented scheme can be 
implemented with current and near future experimental 
techniques. 
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H,s\0)^^u\0)^ + v\l)^ , 

i?eff|l)i=^|l)i+HO)i , 



(A8) 
(A9) 



Solving the corresponding equations for u and v one finds 



u = f2, 
' = -7r92 : 



(AlO) 
(All) 



APPENDIX A: ADIABATIC ELIMINATION FOR 
THE Xl-GATE 

The derivation of Hamiltonian Hx , Eq. (|25|l , is given 
in more detail. The starting point is the Hamiltonian 
Eq. IjlSI) for three atoms in three lattice sites. Since the 
number of atoms in state a and b is conserved the states 
decouple in two subspaces where the first one contains 
only states with one atom in state a and two atoms in 
state b and the other one contains two atoms in state a 
and one atom in state b. The states with three atoms in 
the same internal state do not occur in our system. For 
simplicity we only consider the second subset, the results 
for the other subset can be derived in the same way. 

Let Q be the projection operator on all states with 
more than one atom in a lattice site, P = 1 — Q. The 
Hamiltonian after adiabatic elimination in first order per- 
turbation theory is given by 

Hc« ^ lim PHP - PHQ{QHQ + el)"igiJF . (Al) 

The inversion of QHQ + el was calculated up to first 
order in the perturbation parameter J/U since in this 
treatment higher orders of J/U are neglected anyway. 
In the basis {|100) , |010) , |001)} the adiabatically elimi- 
nated Hamiltonian has the form 
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Ua {{Jt^ 


)2 + (j(^')2)+4C/„,(j1"))2 



UaUab 



(A2) 

(A3) 
(A4) 

(A5) 



/i - /2 + 251 , 

92 = 251 , 
/a = /2 + ffi • 



(A12) 
(A13) 
(AM) 



Because of the symmetry between the states with two 
atoms in state a and one in state b and vice versa the only 
physical solution for these equations is given by Eqs. ^ 



APPENDIX B: QUBITS REMAIN IN THE DFS 

In this appendix we show that the DFS qubits remain 
decoherence free in the case without any interaction be- 
tween the physical qubits, which is relevant for the fast 
i?^(7r) rotation and the CPHASE gate. Since in this case 
there can be more than one atom in a lattice site we in- 
troduce Fock states |n^, n^; n^, n^; n^, n^), where 

(n^) is the number of atoms in mode a (b) in lattice site 
j — 1, 5 , respectively. 

The fast Rl{n) 

ting Ua,b = Uab 



rotation can be performed by set- 



0, J. 



(a.b) j{a.h) j{a,b) 



J 



0. 



In this case the corresponding time evolution operator 
U{t) — exp{—iHt/h), where H is given by Eq. (|15|l . can 
be calculated analytically and we get 



uit) |0)^ = |0)^ , 



(Bl) 



for all t and thus the system never leaves the DFS. The 



state |1)^ evolves as 
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U{t) |1)^ = cos ( ^) |1)^ + -1= sin { \ [|0, 2; 0, 0) + |0, 0; 0, 2)] ® |1, 0; 0, 0; 1, 0) 



h J ' 2%/3 V ^ / I \/2 

+ -^[|2,0;0,0) + |0,0;2,0)]® |0,1;0,0;0,1)- [|1,1;0,0) + |0,0;1,1)]® [|l,0;0,0;0,l) + |0,l;0,0;l,0)]j 

(B2) 



Clearly, the first term cos{2Jt/h) |1)^ remains in the DFS. In order to see that this is also true for the second part we 
rewrite the states in first quantization. Let |q;/3)^ denote a state where the first atom is in state |a) and the second 
one in state |/3), independent in which lattice site they are stored. Due to symmetrization we get 



-i=[|l,l;0,0) + |0,0;l,l)]^i=[|01)^ + |10)J , (B3) 
-i=[|2,0;0,0) + |0,0;2,0)]^|00)^, (B4) 
-i=[|0,2;0,0) + |0,0;0,2)]^|ll)^. (B5) 



By rewriting Eq. ljB2|l we get 



U{t) |1)^ = cos ( |1)^ + sin (^^) {2 |11)^ ® 100)^ + 2 |00)^ ® |11)^ 



ft ; ' 2V3 V ft / I /A I /A ^gg^ 

[|01)^ + |10)J®[|01)^ + |10)J} . 



The part proportional to the sine function is again decoherence free, because it is just the state ll)^^. 

The case of the CPHASE gate is more complicated because the atoms can hop into an additional lattice site. The 
time evolution of state |0)^ is given by 



Uit) |0)^ - [|1, 0; 0, 1) - |0, 1; 1, 0)] ® jcos \y—f-j «; 0, 0; 0, 1) - |0, 1; 0, 0; 1, 0)] 



^Vll, 

(B7) 



-^sin|^^^ [|1,0;0,1;0,0) -|0,1;1,0;0,0) + |0,0;1,0;0,1)-|0,0;0,1;1, 
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For we get 

U{t) |1)^ = [|1, 0; 1, 0) + |0, 1; 0, 1)] ® | ^ + cos (^^^ j j [|0, 1; 0, 0; 1, 0) + |1, 0; 0, 0; 0, 1)] 

+ ^ sin (^^^ j [|0, 0; 0, 1; 1, 0) + |0, 0; 1, 0; 0, 1) + |0, 1; 1, 0; 0, 0) + |1, 0; 0, 1; 0, 0)] 
- ^sin^ (^^j [|1,1;0,0;0,0)+2|0,0;1,1;0,0) + |0,0;0,0;1,1)]| 
+ |O,1;O,1)0|^ |^3 + cos|^^^jj |1,0;0,0;1,0) 

+ ^sin|'?^'j [|1,0;1,0;0,0) + |0,0;1,0;1,0)] (B8) 



1 



: sin^ I I [|2, 0; 0, 0; 0, 0) + 2 |0, 0; 2, 0; 0, 0) + |0, 0; 0, 0; 2, 0)] 



2^2 V ^ / 
+ |1,0;1,0)® |i (^3 + cos|^^^^^ |0,1;0,0;0,1) 

+ ^ sin 1^^^ j [|0, 1; 0, 1; 0, 0) + |0, 0; 0, 1; 0, 1)] 

1 ., sin^ ( ^ 1 [|0, 2; 0, 0; 0, 0) + 2 |0, 0; 0, 2; 0, 0) + |0, 0; 0, 0; 0, 2)] 



2V2 V ^ / 



Making analogous substitutions as above we see that thes e qubits are again decoherence-free at all times. 
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